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ABSTRACT 

We conduct hydrodynamical numerical simulations and find that precessing 
massive slow jets can inflate fat bubbles, i.e., more or less spherical bubbles, that 
are attached to the center of clusters of galaxies. To inflate a fat bubble the 
jet should precess fast. The precessing angle 9 should be large, or change over 
a large range < 9 < # max ^ 30 — 70° (depending also on other parameters), 
where 9 = is the symmetry axis. The constraints on the velocity and mass 
outflow rate are similar to those on wide jets to inflate fat bubbles. The velocity 
should be vj ~ 10 4 km s" 1 , and the mass loss rate of the two jets should be 
2Mj ~ 1 — 5OM yr _1 . These results, and our results from a previous paper 
dealing with slow wide jets, support the claim that a large fraction of the feedback 
heating in cooling flow clusters and in the processes of galaxy formation is done 
by slow massive jets. 



1. INTRODUCTION 

Many of the X-ray deficient bubbles in galaxies and clusters of galaxies reside very close 
to the center of the cluster (or galaxy) and are fully or partially surrounded by a dense shell, 
e.g., Perseus (Fabian et al. 2000), Abell 2052, (Blanton et al. 2003), Abell 4059 (Heinz et 
al. 2002), and HCG 62 (Vrtilek et al. 2002; Morita et al. 2006). We term these more or less 
spherical bubbles 'fat bubbles'. Fat bubbles are defined by the following properties: (1) Fat 
bubbles come in pairs, each on opposite sides of the equatorial plane. In some cases there is 
departure from axisymmetry, and they are not exactly opposite. (2) They touch each other 
at the center, and by that form an hourglass structure (like the figure '8'). One bubble of 
the hourglass structure is referred to as a fat bubble. (3) The density inside the bubbles is 
much lower than that of their surroundings (ambient gas). (4) They are fully or partially 
surrounded by relatively thin shell that is denser than the surroundings. (5) In some cases 
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their boundary, on the far side from the center, is open. In classifying planetary nebulae, for 
example, this structure is refereed to as a bipolar nebula. There are many such planetary 
nebulae that are well resolved in the visible band, and the bipolar structure is well defined. 
In many cases there are similar structures in planetary nebulae and in clusters of galaxies 
(Soker & Bisker 2005). The best examples of the hourglass type of structure we aim to study 
are the bubbles in Perseus (Fabian et al. 2000) and in A 2052, (Blanton et al. 2001). Other 
cases are cited above and in Soker & Bisker 2006. If the bubble rises through the ICM and 
do not touch each other any mode, then they are not 'fat' any more because they don't form 
an hourglass structure anymore. We estimate that ~ 30 — 50% of the bubbles are, or were, 
fat bubbles. 

In recent years, two and three-dimensional hydrodynamical simulations of jets and bub- 
bles in clusters of galaxies, were conducted to study different aspects of their interaction with 
the intra-cluster medium (ICM), such as heating the ICM (e.g., Basson & Alexander, 2003; 
Heinz & Churazov, 2005; Reynolds et al., 2005; Heinz et al., 2006; Vernaleo & Reynolds, 
2006; Binney et al., 2007; Ruszkowski et al., 2007; Alouani Bibi et al., 2007; Briiggen et al., 
2007). In our studies we aim to understand the conditions leading to the formation of fat 
bubbles. 

In previous papers (Soker 2004, 2006; Sternberg et al. 2007, hereafter Paper I) we 
proposed that in order to inflate fat bubbles, either the jet's opening angle has to be large, 
i.e., wide jets, or the jet is narrow but its axis has to change its direction. The change in 
direction can result from precession (Soker 2004, 2006), random change (Heinz et al. 2006), 
or a relative motion between the ICM and the active galactic nucleus (AGN; Loken et al. 
1995; Soker & Bisker 2006; Rodriguez-Martmez et al. 2006). 

In Paper I we conducted two-dimensional hydrodynamical simulations of wide jets ex- 
panding in the ICM. We found that wide jets can indeed inflate fat bubbles that reside very 
close to the center of the cluster. For this to occur, we found that the jets should have high 
momentum flux. Typically, the half opening angle should be a > 50°, and the large momen- 
tum flux requires a jet speed of vj ~ 10 4 km s _1 , i.e., highly non-relativist ic, but supersonic 
(vj c± 3000 — 3 x 10 4 km s _1 ). Narrow relativistic jets can exist in parallel with the slow wide 
outflow, but they will not lead to the inflation of fat bubbles. The inflation process involves 
large vortices and local instabilities which mix some ICM with the hot bubble. These results 
predict that most of the gas inside the bubble has a temperature of 3 x 10 8 < T5 < 3 x 10 9 K, 
and that large quantities of the cooling gas in cooling flow clusters are expelled back to the 
intra-cluster medium and heated up (Soker & Pizzolato 2005). In Paper I we suggested that 
the magnetic fields and relativistic electrons that produce the synchrotron radio emission 
might be formed in the shock wave of the non-relativistic jet. Motivated by our earlier 



-3- 



results, in this paper we examine the inflation of fat bubbles by narrow precessing jets. 



The simulations were performed using the Virginia Hydrodynamics- 1 code (VH-1; Blondin 
et al. 1990; Stevens et al. 1992), as described in Paper I. We simulated a three-dimensional 
axisymmetric flow, so practically we simulated a quarter of the meridional plane using the 
two-dimensional version of the code in spherical coordinates. The symmetry axis of all plots 
shown in this paper is along the horizontal axis (the x axis) while the equitorial plane is ver- 
tical. Radiative cooling and gravity were not included, since the total time of the simulation, 
tsim ~ 10 7 yr, is somewhat shorter than the gravitational time scale, and much shorter than 
the radiative cooling time. This preliminary report aims to emphasize the jet properties that 
determine whether or not the required bubble is inflated, hence, these omissions are justified. 

We used the (3 model (with (3 = 1/2) as the initial density profile of the ICM, 



with p c = 2.16 x 10 -25 g cm -3 and r = 100 kpc (see Paper I and references therein). The 
ICM temperature is 2.7 x 10 7 K. The box size used in our simulations was 30 kpc x 30 kpc 
(one quarter of the meridional plane). We used a 128 x 128 evenly spaced grid. As discussed 
in Paper I, higher resolution does not change the large scale behavior, thus, this resolution is 
sufficient for our study. In Figure [1] we present density maps (logarithem of the density. The 
density is given in g cm -3 .) of two runs with the same parameters but different resolution 
(128 x 128 and 256 x 256). This figure clearly shows that the large scale behavior does 
not change significantly with the increase in resolution. Moreover, we note that in the 
zone relevant to our simulations, r < 30 kpc, the assumption of a constant temperature is 
reasonable (see fig 1. of Pizzolato & Soker 2005). 

The narrow jet was injected at a radius of 0.1 kpc, with constant mass flux Mj (per one 
jet) and a constant radial velocity inside a half opening angle a = 5°. Therefore, the total 
kinetic power of one jet is Ej = MjVj/2. . The symmetry axis of the jet is at an angle 9(t) in 
respect to the symmetry axis of the problem (the x axis) . This is the precessing angle which 
is a function of time. Because of the axisymmetric nature of our problem, the meaning of a 
precessing jet in these simulations is that the narrow jet precess around the symmetry axis 
very rapidly. Namely, the precessing period around the symmetry axis is much shorter than 
any other relevant time scale in the problem, e.g., the time scale over which 9 is changing. 

Due to the numerical nature of the jets injection, in some cases a fraction of the injected 
mass does not succeed in mounting our grid and remains in the first cells of the grid. In the 
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next time step the jet properties are reinserted into these cells and this mass is lost. Moreover, 
if the jet is precessing, then the movement of the injection zone my increase/decrease the 
effective opening angle of the jet, in effect changing Mj, such that only after the simulation 
we can determine the mass loss rate and kinetic luminosity of the jet. The kinetic luminosity 
of the jet is constant, and is equal to 5PV/2t, where P is the pressure and V is the volume 
inflated after time t. We calculated this value at several times for each case, and was taken to 
be the kinetic luminosity. All tables and plot captions specify the estimated jet luminosity. 
We consider three basic types of precessing jets: 

(i) Fixed precessing angle, i.e., 9 is constant. 

(ii) A constant rate of change in the precession angle, i.e., at constant time interval d6/dt 

is constant for 6 < # max . 

(iii) Random precession, i.e., the jet axis has the same probability to take any direction 
within a maximum angle # max - This is done by taking d(cos9)/dt to be constant. 

We studied 12 cases per each precession type (fixed, constant rate change, and random). 
The parameters of the different cases are given in Tables [IE 

We used slow massive jets, as have been used before in a number of numerical studies 
(e.g., Paper I; Alouani Bibi et al. 2007). We further discuss the usage of slow massive jets 
in section 6. 



3. RESULTS: PRECESSION AT A CONSTANT ANGLE 

In Figure [2] we show the density map (logarithem of the density. The density is given 
in g cm -3 .) at different times for a jet with a fixed precession angle of 6 = 15° in respect 
to the symmetry axis (Fix2, see Table []]). The jet has a half opening angle of a = 5°, 
Vj = 7750 km s -1 , and Ej ~ 1.1 x 10 44 erg s -1 . For this case the mass injection rate into 
one jet is Mj ~ 6M yr -1 . Namely, the two jets expel mass back to the ICM at a high 
rate of ~ 12M yr -1 . At first the jet inflates a cavity of low density matter in the shape of 
a torus. At later times the jet is bent towards the symmetry axis, thereafter, it continues 
it's propagation in a manner quit similar to that of a non precessing jet with a ~ 20° half 
opening angle (see Figure 1 model 2 in Paper I). In this case, it is obvious to see that a fat 
bubble was not inflated. The arrows in the plot represent the velocity of the flow. For visual 
clarity we divided the velocities into groups, each represented by an arrow of a predetermined 
length: 
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Fig. 1. — Density maps for two runs with the same parameters but different resolution 
(128 x 128 and 256 x 256) shown at the same time (t = 27 Myr). The small scale behavior 
might differ, but the large scale behavior, which is what we are interested in, does not differ. 
We emphasize that the x axis is the symmetry axis. In all figures the density is given in 
g cm -3 and in log scale. 



Precession at a fixed angle 


Run 


O(deg) 






Morphology 


Fig 


ure 


Fixl 


15° 


7750 


0.24 


Thin and narrow jet and cocoon 






Fix2 


15° 


7750 


1.1 


Thin and narrow jet and extensive cocoon 


Fi { 


m 


Fix3 


15° 


23250 


0.15 


Very Thin and narrow jet and cocoon 






Fix4 


15° 


23250 


1.8 


Thin and narrow jet and extensive cocoon 


Fi { 




Fix5 


30° 


7750 


0.14 


Fat Bubble 


Fi { 




Fix6 


30° 


7750 


1.8 


Fat bubble 


Fi { 




Fix7 


30° 


23250 


1.2 


First torus then spherical cavity at center of cluster 


Fi { 




Fix8 


30° 


69750 


1.2 


First torus then an ellipsoid cavity at center of cluster 






Fix9 


60° 


7750 


0.2 


First torus then an ellipsoid cavity at center of cluster 






FixlO 


60° 


7750 


1.8 


Both tori merge to create a 'doughnut' 


Fi { 




Fixll 


60° 


23250 


0.14 


First torus then an ellipsoid cavity at center of cluster 






Fixl2 


60° 


23250 


2.2 


First torus then spherical cavity at center of cluster 







Table 1: Parameters of the fixed angle precession runs, where, 9 is the angle (in degrees) 
between the symmetry axis of the jet and the symmetry axis of the problem, Vj is the jet 
velocity in Km s _1 , Lj is the jet luminosity in erg s _1 (per one jet). The cluster sound speed 
is 775 Km s _1 . 
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Constant rate of change in the precession angle 


Run 


T(Myr) 






Morphology 


Figure 


Conl 


0.1 


7750 


1.9 


Fat bubble 


Fig El 


Con2 


0.1 


23250 


1.1 


Thin jet shedding large vortecies 


Fig El 


Con3 


0.1 


69750 


1.8 


Very thin jet and extensive cocoon 


- 


Con4 


1 


7750 


2 


Narrow and elongated cavity (not a fat bubble!!!) 




Con5 


1 


23250 


1.5 


Narrow and elongated cavity 




Con6 


1 


69750 


1.2 


Narrow and elongated cavity 




Con7 


5 


7750 


1.6 


Narrow and elongated cavity extensive backflow 


Fig El 


Con8 


5 


23250 


1.6 


Narrow jet with extensive backflow 




Con9 


5 


69750 


2 


Narrow jet with extensive backflow 




ConlO 


30 


7750 


1.8 


Elongated clumpy cavity, extensive backflow 




Conll 


30 


23250 


1.5 


Narrow jet with extensive clumpy cocoon 




Conl2 


30 


69750 


1.6 


Narrow jet with narrow clumpy cocoon 





Table 2: Parameters of the constant rate of change in precession angle runs. Same parameters 
as in Table [U In addition, T is the precession period in Myr. 



Random change in the precession angle 


Run 


T(Myr) 




L,(10 44 ^) 


Morphology 




^ure 


Rani 


0.1 


7750 


1.4 


Fat bubble 


Fi| 




Ran2 


0.1 


23250 


2 


Thin jet shedding very large vortecies 






Ran3 


0.1 


69750 


1.5 


Thin jet and cocoon, extensive backflow 






Ran4 


1 


7750 


1.8 


Elongated cavity 






Ran5 


1 


23250 


1.6 


Elongated cavity 






Ran6 


1 


69750 


1.5 


Narrow and elongated cavity extensive backflow 






Ran7 


5 


7750 


1.5 


Elongated clumpy cavity 






Ran8 


5 


23250 


1.8 


Elongated clumpy cavity 






Ran9 


5 


69750 


2 


Elongated clumpy cavity, extensive backflow 






RanlO 


30 


7750 


1.8 


First torus then fat bubble 


Fk; 




Ranll 


30 


23250 


1.5 


First torus then elongate cavity, extensive backflow 






Ranl2 


30 


69750 


1.6 


Torus then narrow jet with clumpy cocoon 







Table 3: Parameters of the random change in precession angle runs. Same as in Table El 
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(i) 0.1c s < vj < 0.5c s - shortest, 

(ii) 0.5c s < ^ < c s: 

(iii) c s < ^ < 5c s , 

(iv) 5c 5 < Vj < 10c s - longest for the M = 10 case, 

(v) 10c s < < 30c 5 - longest for the M = 30 case, 

where c s = 775 km s _1 is the speed of sound. For visual clarity we also omitted velocities of 
Vj < 0.1c s . This division is true for all plots shown in this paper. 

In Figure [3] we show the density map at different times for a jet at a fixed precession 
angle of 9 = 30° (Fix5, see Tabled]), with a half opening angle of a = 5°, vj = 7750 km s _1 , 
and Ej ~ 1.4 x 10 43 erg s _1 . In a similar manner to the former case, the jet initially inflates 
a torus shaped cavity, and later it is bent towards the symmetry axis. In contrast to the 
former case we see that at t = 45 Myr the jet has inflated a bubble that together with 
the contra-bubble will form a bipolar structure. The bubble in this case is much closer to 
being spherical than bubbles in runs that we consider unsuccessful in forming fat bubbles 
(they don't form a bipolar structure). Even at t = 60 Myr the cavity can still be termed an 
elongated bubble. In contrast to the bubbles we showed in Paper I, there is no flow of low 
density matter towards the equatorial plane, which is in good agreement with observations. 
For the same parameters, but with a luminosity of <~ 2 x 10 44 erg s _1 (Fix6, see Table [1]), we 
also got a fat bubble of low density, as shown in Figure 0J In contrast to the lower luminosity 
case, we see a non-negligible flow of low density matter towards the equatorial plane, which 
is not in very good agreement with observations. At early time a torus, rather than a fat 
bubble, is inflated. At present, we know of no observations of such bubbles. 

In most of the cases that we ran we ended up with either a narrow jet, as shown in 
upper panel of Figure [5] (Fix4,see Table [1]), or with a fat bubble with a significant flow of 
low density matter to the equatorial plane (in effect a spherical or ellipsoid cavity at the 
center of the cluster), as shown in the middle panel of Figure [5] (Fix7, see Table [1]). In some 
of the cases with 9 = 60° (Fixl0,see Table [1]) the torus inflated by the simulated jet merged 
at the equatorial plane with the torus of the unsimulated jet (we remind the reader that we 
imposed reflecting boundary conditions in the equatorial plane) . In effect we got a large low 
density torus, with an elliptic cross section, at the equator as shown in the lower panel of 
Figure For the sake of clarity we state that all of these cases are not in agreement with 
observations. Therefore, large constant precessing angle 9 > 40°, cannot form the observed 
fat bubbles. In all cases the typical temperature of the low density gas in the cavities or 
cocoon was 10 8 K < T b < 10 9 K. 
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Fig. 2. — Density maps (logarithem of the density [g cm -3 ]) for a jet at a fixed precession 
angle 6 = 15° (Fix2, see Tabled]), given at three different times (t = 2.5, 5, and 20 Myr). The 
jet has a half opening angle of a = 5°, an injected velocity of vj = 7750 km s -1 , and power 
of (one jet) £j ~ 1.1 x 10 44 erg s -1 . Only one quarter of the meridional plane is showed, 
as the other three are symmetric to it. The x-axis (horizontal) is the symmetry axis, while 
the y-axis (vertical) is in the equatorial plane. The arrows represent the velocity of the flow: 
0.1c s < vj < 0.5c 5 (shortest), 0.5c s < vj < c s , c s < vj < 5c s , and 5c 5 < vj < 10c s (longest in 
this case). 
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Fig. 3. — Density maps for a jet at a fixed precession angle of 9 
with a half opening angle of a = 5°, vj = 7750 km s" 



30° (Fix5, see Tabled]), 
and Ej ^ 1.4 x 10 43 erg s" 1 . The 
cavity inflated by the jet reaches a more or less spherical shape at t ~ 45 Myr, and remains 
spherical, though a little elongated even at t = 60 Myr. The equatorial plane (vertical in the 
figures) remains devoid of low density matter, in good agreement with observations. There 
are no numerical fluctuations in the equatorial boundary either, unlike along the symmetry 
axis (horizontal axis) The arrows represent the velocity of the flow as in Figure [3 
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4. RESULTS: A CONSTANT RATE OF CHANGE IN THE PRECESSION 

ANGLE 

Figure [6] shows the density maps for a precessing jet with a constant rate of change in 
the precession angle (i.e., 9 oc £, where the symmetry axis of the jet changes in the range 
5° < 9 < 65°). We remind the reader that in all cases presented here, the jet is assumed to 
precesses rapidly around the symmetry axis, so it rotates in the (j) direction many times while 
9 is being changed; the (j) coordinate is not calculated in the simulations. The case shown 
in this figure has a precession period T prec = 0.1 Myr (Conl, see Table [2]). This period is 
much shorter than the typical expansion time of the bubble formed. As a result of that the 
narrow jet's interaction with the ICM resembles that of a wide opening angle jet with a half 
opening angle of ~ 70° (Paper I). For comparison, in figure [3 we show the case of a wide 
angle jet taken from Paper I. As can be seen in figure El at short times of (t < 5 Myr) the 
low density bubble is more or less spherical and there is almost no flow of low density matter 
to the equatorial plane. At longer times of t ~ 10 Myr, the flow of low density matter to the 
equatorial plane is substantial, though the cavity itself can still be termed a fat bubble. We 
conclude that for these parameters a fat bubble is formed. The difference in the volume of 
the bubbles is due to the fact that the actual M of the wide angle case was slightly smaller 
then that of the precessing case, and therefore the bubble inflated by it was slightly smaller 
(this is due to the numerics associated with the jet injection, as elaborated in section EJ). 

Figure [8] shows the case with higher jet velocity of vj = 23250 km s _1 , T = 0.1 Myr, 
Ej ~ 1.2 x 10 44 erg s _1 (Con2,see Table [2]). As with wide jets (Paper I), fast jets of Vj > 
20, 000 km s _1 (the exact speed limit depends on the other parameters) do not form fat 
bubbles, but rather the jets propagated close to the axis of symmetry (the x axis) in the 
manner of a narrow jet with an extensive cocoon. 

In the case of constant rate of change in the precession angle, the jet spends equal time 
in small and large precession angles. At small precession angle the jet has more momentum 
per unit area, making it easier for the jet to break through the denser ICM. This can be seen 
in simulation with longer precession periods, e.g., T prec > 1 Myr (not shown here) where 
at small precession angles the jet's propagation is easier than at larger angles, resulting in 
a propagation along the symmetry axis, i.e., in the resemblance of a narrow jet, as can be 
seen for example in Figure [9] (Con7, see Table [2]). No fat bubble is formed for a too long 
precession period. In all cases simulated in this section the typical temperature of the low 
density gas in the cavities or cocoon was of the order 10 8 K < T5 < 10 9 K. 
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5. RESULTS: RANDOM PRECESSION 

A more physically acceptable change in the precession angle is one in which the jet covers 
a constant solid-angle per unit time, i.e., d(cos9)/dt is constant. Because the precession 
period about the symmetry axis is assumed to be very short, the jet spends longer periods 
of time at larger angels, in respect to the symmetry axis, than at small ones. The shorter 
time spent at small angels reduces the break-through period along the symmetry axis , as 
experienced in the uniform change in precession angel case (section [3]) . The reduction in the 
break-through period allows more matter to spread farther from the symmetry axis, in effect 
inflating a fat bubble of low density matter. This can be seen in Figure [10] which shows 
the density map of a randomly precessing jet with a precession period of T prec = 30 Myr 
(RanlO, see Table [3]). The other parameters are vj = 7750 km s _1 , Ej ~ 1.8 x 10 44 erg s _1 , 
and the boundary of the precession angle (of the jet's axis) are 5° < 9 < 45°. The jet starts 
its precession off- axis and inflates a toroidal cavity as seen at t = 10 Myr. At t = 20 Myr 
we see that the cavity is in the shape of a fat ellipsoid bubble. At t = 25 Myr the bubble is 
more or less spherical with a radius R ~ 20 kpc. There is a flow of low density matter to 
the equatorial plane, but it is minimal. 

As in the constant rate of change case, the random precessing jet with a short precessing 
period, T prec = 0.1 Myr (Rani, see Table BJ) resulted in an interaction with the ICM similar 
to that of a wide jet with a half opening angle of a ~ 50°. Figure [TT] shows the density maps 
for this case. All along the simulation the cavity inflated by the jet is more or less spherical, 
though at t = 25 Myr it gets a little elongated. In contrast to observations there is a flow of 
low density matter to the equatorial plane. 

For higher velocity jets (vj > 2 x 10 4 km s _1 ) the results (not shown here) of the 
interaction between the jet and ICM are the propagation of the jet close to the symmetry 
axis (i.e., a narrow jet with a cross section radius of 1 — 2 kpc) with an extensive cocoon of 
low density matter shed by the propagating jet in the form of vortices. The cross section 
radius, of the jet and the cocoon, is typically 3 — 5 kpc. 

Over all, we conclude that rapidly and randomly precessing massive slow jets can inflate 
fat bubbles, similar to those inflated by wide jets (Paper I). In all cases the typical temper- 
ature of the low density gas in the cavities or cocoon was of the order 10 8 K < T5 < 10 9 K. 

6. DISCUSSION AND SUMMARY 

We showed that precessing slow jets can inflate fat bubbles attached to the center of 
the galaxy clusters. By slow jets we refer to supersonic but highly non-relativistic jets. In 
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our axisymmetrical simulations the 3D problem can be simulated with a 2D grid, and the 
simulations were of jets precessing rapidly in the (j) coordinate around the symmetry axis 
(the (j) coordinate is not included in our simulations). Namely, the rotation time of the jet's 
axis around the symmetry axis (the x axis in our figures) is much shorter than the time over 
which the precessing angle 9 varies, the length of our simulations was chosen to match the 
10 — 50 Myr age of most observed bubbles (Birzan et al. 2004; McNamara & Nulsen 2007). 

The main criteria we find for the inflation of fat low density bubbles by precessing jets 

are: 

1. The jet's velocity should be vj ~ 10 4 km s _1 . Using our results obtaned here and our 
results from Paper I, we conclude that the range over which slow jets can inflate fat 
bubbles is 3000 km s _1 < v 3 < 2 x 10 4 km s _1 . Jets with higher velocities form narrow 
expanding jets with extensive cocoons. These jet velocities form bubbles with interior 
temperatures of 10 8 K < T b < 10 10 K. 

2. For a jet with a fixed precession angle 9 (measured from the the symmetry axis; x in our 
figures), this angle should not be too small or too large. A small angle causes the jet to 
propagate along the symmetry axis, and a large angle leads to the bending of the jet 
towards the equatorial plane. The constraint on the precessing angle is 30° < 9 < 50° 

3. The jet should spend more time at large angles in order to reduce the break-through 
along the symmetry axis that happens when it's precessing angle is small. This favors 
a case with random precession, i.e., the jet axis covers constant solid angle per unit 
time and the precessing angle is bound in the range 0° < 9 < # max ~ 50°. 

4. For any prescribed precession behavior and parameters that can form a fat bubble, the 
maximum precessing angle should be quite large, 9 max ^ 30 — 70°. 

5. The slow velocity and large energy of the jet inflating fat bubbles require that the two 
opposite jets carry large amount of mass. The two jets together can expel back to the 
ICM a mass at a rate of M back ~ 1 — 5OM yr" 1 . 

Let us elaborate on these points, and on the physics behind them. As discussed in 
previous papers (Soker 2004; 2006; Paper I), the basic condition for a jet to inflate a fat 
bubble is that the jet's head will reside inside the bubble, or, if the jet's head is outside the 
bubble, that the jet's head will not "run away" from the expanding bubble. For example, in 
the case of constant large precessing angle (Run 'FixlO' in Figure Ej), the jet's head revolves 
around the symmetry axis but at a large distance. It inflates a local low density region, a 
torus around the axis. But this torus expands too slowly to inflate a fat bubble. Namely, 
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the motion of the jet's head is too fast for the expanding shocked gas. This condition for 
the jet's direction not to escape from the expanding bubble formed by its shocked material 
is given in equation (17) of Soker (2006). This explains why the precessing angle cannot be 
too large (point 2 above). 

The same principle holds for the jet's head not to expand too fast in the radial direction. 
If the expansion of the jet is concentrated within a small solid angle, e.g., a small precession 
angle, the large momentum flux (ram pressure) will result in a jet's head that moves radially 
faster than the expansion of the shocked gas. A long bubble will be formed instead of a fat 
bubble. This explains points 3 and 4 above and the lower limit on the precessing constant 
angle in point 2. Qualitatively, this condition on the solid angle of the expanding jet is given 
by equation (14) in Soker (2004) that was derived for a wide jet instead of a precessing jet, 
and was shown to hold in the numerical simulations of wide jets (Paper I). 

Regarding point 1 above. If the jet's speed is too low, for a given energy it is too dense. 
As a result of that, the bubble expands too slowly and the jet's head moves radially too 
fast. This is also seen by equation (14) in Soker (2004), and was shown to hold for wide jets 
(Paper I). If the jet's material is too fast, the shocked jet's gas expands much faster than 
the jet's head. In that case the shocked low density gas fills the entire inner region, and 
the two opposite jets form a large elliptically-shaped bubble instead of a bipolar (hourglass) 
structure of two fat bubbles. This is very similar to the case of wide angle (Paper I). Note 
that in some of the simulations no dense ICM gas is left in the equatorial plane near the 
center (e.g., Figure El see also Paper I). This shows that dense ICM in the equatorial plane 
near the center we find in some runs (e.g., Figures ITOl see also Paper I) is real and not a 
numerical artifact. 

Point 5 has far reaching implications. Are such a high mass outflow rates as we find 
here and in Paper I where we simulated wide jets, Mback — 1 — 5OM yr _1 , compatible 
with observations? Are wide jets simulated in Paper I (or wide precession angle simulated 
here) compatible with observations? The answer to both questions seems to be positive. In 
Paper I we have already discussed indications for AGNs that blow slow jets, some of them 
with wide angles (Crenshaw & Kraemer 2007; Behar et al. 2003; Kaspi & Behar 2006). 
These recent observational results (see also de Kool et al. 2001), and our numerical results 
support the model where the feedback in both cooling flow clusters and in the process of 
galaxy formation occurs mainly (but not solely, as relativistic narrow jets also exist) by slow 
massive jets, as suggested and discussed by Soker & Pizzolato (2005). The massive jets 
imply that not only energy, but mass as well is part of the feedback cycle (Soker & Pizzolato 
2005; Pizzolato & Soker 2005). Massive jets were also considered before by, e.g., Begelman 
& Celotti (2004) and Binney (2004), and were simulated by Omma et al. (2004), who took 
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the jet speed and mass outflow rate to be Vj = 10 4 km s" 1 and 2M yr -1 , respectively. In 
the simulation of Heinz et al. (2006) one jet has a mass loss rate of 35M yr _1 . This implies 
that the two jets inject 7OM yr _1 into the ICM. Without stating it, Heinz et al. (2006) 
followed the suggestion of Soker & Pizzolato (2005), that a large fraction or even most of 
the gas that cools to low temperatures in cooling flow clusters gains energy directly from the 
central black hole, and is injected back to the ICM. 

The main effect discussed in the present paper is a dynamical effect with supersonic 
velocities. Therefore, the temperature (and entropy) profile of the ICM has little effect on 
the conclusions. In a forthcoming paper we will follow the bubbles as they buoy to larger 
distances (~ 100 kpc). There the entropy will have a crucial role, and a more realistic 
temperature profile will be used. 

There are some discrepancies between our results and observations. 

(1) In many of the simulated cases, but not in all of them, a non negligible low density 

matter flows to the equatorial plane and stays there. This is generally not observed. 
It is quite possible that magnetic tension will suppress this low density backflow. This 
is because the backflow is expected to stretch magnetic field lines. In addition, gravity 
might cause this low density equatorial gas to buoy outward. Gravity and magnetic 
fields are not included in our simulations. We will include gravity in future simulations. 
However, neither gravity nor magnetic fields are expected to prevent the presence of 
hot low density gas in the equatorial plane in cases of very fast jets where the backflow 
toward the equator is strong. 

(2) In some cases a torus shaped cavity is inflated at the beginning of the simulation, 

achieving a spherical shaped cavity only at later times. We remind the reader that we 
assumed a very short-period precession about the symmetry axis, in effect reducing 
the simulations dimensions to two. A full three dimensional simulation of a randomly 
precessing jets (and even a pulsed jet, whose activity is turned on and off), might result 
in the inflation of low density small bubbles and not a torus as in the 2D simulations. 
These small bubbles merge on a short time to form a larger bubble. We leave the 
investigation of this process to a full 3D simulations in the future. 

Concerning these discrepancies and the observations discussed above, we suggest that 
in most, but probably not in all, cases fat bubbles are inflated by wide jets or rapidly and 
randomly precessing jets. 

We thank John Blondin for his immense help with the numerical code. We thank Ehud 
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This research was supported by the Asher Fund for Space Research at the Technion. 
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Fig. 4. — Density maps for a jet at a fixed precession angle of 9 = 30°, with a half opening 
angle of a = 5°, v 5 = 7750 km s" 1 (Fix6, see Table OQ), and Ej ~ 2 x 10 44 erg s" 1 . The 
jet inflates a more or less spherical cavity. The bubble is shown at t = 40 Myr. There is a 
non-negligible flow of low density matter towards the equatorial plane, which is not in good 
agreement with observations. The arrows represent the velocity of the flow as in Figure [3. 
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Fig. 5. — Density maps for three different cases of jets at a fixed precession angles with a 
half opening angle of a = 5°. In the upper panel (case A: Fix4,see Tabled]) we show a high 
velocity jet (vj = 23250 km s _1 ) with a small precession angel (9 = 15°), and jet luminosity 
of Ej c± 1.8 x 10 44 erg s _1 . In the middle panel (case B: Fix7,see Table HD we show high 
velocity jet (vj = 23250 km s _1 ) with an intermediate precession angel (9 = 30°), and jet 
luminosity of Ej ~ 1.2 x 10 44 erg s _1 . In the lower panel (case C: Fixl0,see Tabled]) we 
show a low velocity jet (vj = 7750 km s _1 ) with a large precession angel (9 = 60°) and jet 
luminosity of Ej ^ 1.8 x 10 44 erg s _1 . In all three cases a fat bubble was not inflated. 
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Fig. 6. — Density maps for a jet with a uniformly changing precession angle, 5° < 9 < 65°, 
with a precession period of T prec = 0.1 Myr, half opening angle of a = 5°, vj = 7750 km s _1 , 
and Ej ~ 1.9 x 10 44 erg s _1 (Conl,see Table [2]). The arrows represent the velocity of the 
flow as in Figure El 
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Fig. 7. — A case of a wide jet (a = 70°) instead of a precessing jet, taken from Paper I. The 
parameters of this run are: vj = 7750 km s _1 , and Ej ~ 0.65 x 10 44 erg s _1 . This plot shows 
the bubble at t = 5 Myr. The arrows represent the velocity of the flow as in Figure El 
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Fig. 8. — Like figure [6] but for a faster jet with Vj = 23250 km s x , i.e., 5° < 9 < 65°, 
T prec = 0.1 Myr, a = 5°, and E 3 ~ 1.2 x 10 44 erg s" 1 (Con2,see Table Ej). The arrows 
represent the velocity of the flow, according to the five velocity ranges given in section 3. 
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Fig. 9. — Density maps for a jet with a constant rate of changing precession angle, 5° < 
9 < 65°, with a precession period of T prec = 5 Myr, half opening angle of a = 5°, vj = 
7750 km s _1 , and Ej c± 1.6 x 10 44 erg s _1 (Con7, see Table Ej). The arrows represent the 
velocity of the flow as in Figure El The simulation was started with a precession angle of 
9(0) = 5°. 
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Fig. 10. — Density maps for a randomly precessing jet with T prec = 30 Myr, a = 5°, vj = 
7750 km s _1 , and Ej ~ 1.8 x 10 44 erg s _1 (RanlO, see Table [3]). The jet's axis precess 
between # min = 5° and 6^ max = 45°. The arrows represent the velocity of the flow as in Figure 
[3 The simulation was started with a precession angle of 9(0) = 45°. 
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Fig. 11. — Density maps for a randomly precessing jet with precession period of T prec = 
0.1 Myr, jet's opening angle of a = 5°, jet's speed vj = 7750 km s _1 , and one jet power of 
Ej ~ 1.4 x 10 44 erg s _1 (Rani, see Table [3]). The jet's interaction with the ICM is similar 
to that of a wide jet with half opening angle of a ~ 50° (Paper I). The arrows represent the 
velocity of the flow as in Figure El 



